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Abstract 

A light-cone QCD-inspired model, with the mass squared operator consisting 
of a harmonic oscillator potential as confinement and a Dirac-delta interac- 
tion, is used to study the 5-wave meson spectra. The two parameters of the 
harmonic potential and quark masses are fixed by masses of p(770), p(1450), 
J/i/j, ip(2S), i^*(892) and B*. We apply a renormalization method to define 
the model, in which the pseudo-scalar ground state mass fixes the renormal- 
ized strength of the Dirac-delta interaction. The model presents an universal 
and satisfactory description of both singlet and triplet states of S- wave mesons 
and the corresponding radial excitations. 

PACS numbers: ll.10.Ef, ll.10.Gh, 12.39.Hg, 12.39.Ki, 14.40-n 

I. INTRODUCTION 

In the effective light-cone QCD theory JT|^] the lowest Fock component of the hadron 
wave function is an eigenfunction of an effective mass squared operator with constituent 
quark degrees of freedom and parameterized in terms of an interaction which contains a 
Coulomb-like potential and a Dirac-delta term. The Fock-state components of the hadron 
light-front wave function can be constructed recursively from the lowest Fock-state com- 
ponent. The interaction in the mass operator comes from an effective one-gluon-exchange 
where the Dirac-delta term corresponds to the hyperfine interaction. The masses of the 
ground state of the pseudo-scalar mesons and in particular the pion structure || were de- 
scribed reasonably, with a small number of free parameters, which is only the canonical 
number plus one — the renormalized strength of the Dirac-delta interaction. 
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The model was extended to include the confining interaction and used to study the 
splitting of the excited pseudo-scalar states from the excited 3 Si vector meson states as a 
function of the ground state pseudo-scalar mass M. In Ref. [ffl, the Coulomb-like and the 
confining interactions were substituted by a harmonic oscillator potential, which allowed 
an analytic formulation. The parameters of the confining interaction in the mass squared 
operator, were fitted to the 3 Si-meson ground state mass and to the slope of the trajectory 
of excited states with the radial quantum number [H. With the renormalized strength of 
the Dirac-delta interaction fixed by pseudo-scalar masses, it was shown that the n-p mass 
splitting, due to the attractive Dirac-delta interaction, is the source of the splitting between 
the masses of the excited states. A reasonable agreement with the data || was obtained. 
Besides, in light-cone framework (mass squared operator appears in the Hamiltonian) , the 
Dirac-delta plus confining harmonic oscillator potential gave a natural explanation of the 
observation of the almost linear relationship between the mass squared of excited states with 
radial quantum number n ||. This reveals some of the physics that are brought by the work 
of Ref. and showed the relation between the 7r and p spectrum, through the pion mass 
scale which defines the renormalization condition of the model. 

In the present paper, this simple model, with the Dirac-delta interaction acting in the 1 5'o 
channel only and harmonic oscillator potential as confinement, is used to investigate the S- 
wave meson spectra from n-p to 77&-T and make predictions for rj t -9 (we do not study rj-u and 
r)'-4>). Instead of using flavor-dependent parameters, parameters are found for the harmonic 
oscillator potential which is valid for not only light mesons but also heavy mesons. In other 
words, the parameters are universal. It is shown that the linear relationship between the 
mass squared of excited states with radial quantum number is still qualitatively valid even 
for heavy mesons like T. The simply model presents reasonable agreement with available 
data and/or with the meson mass spectra given by Godfrey and Isgur [[7[]. 

This paper is organized as follows. In section [0], we give very briefly the extension of the 
light-cone QCD-inspired theory for which the mass squared operator of a constituent quark- 
antiquark system includes a confining interaction Q. The renormalization of the theory 
using the subtracted equations for the transition matrix || of the model can be found in 
Ref. |3|,|]] thus is omitted here. In the same section, we present the Dirac-delta term plus 
harmonic oscillator potential approach and solve it with the T matrix method developed in 
Ref. [[],f|| . The results and discussion are presented in section [II]]. A brief summary is given 
in the last section. 



II. EXTENDED LIGHT-CONE QCD-INSPIRED THEORY WITH DIRAC-DELTA 
AND HARMONIC OSCILLATOR CONFINING POTENTIAL 

In this section we review our previous work in which we have extended the renor- 
malized effective QCD-theory of Ref. || to include confinement. In the effective theory the 
bare mass operator equation for the lowest Light-Front Fock-state component of a bound 
system of a constituent quark and antiquark of masses mi and 777,2, is described as MM 



A- 



+ m\ 

x 



+ 



1 — X 



2 



dx'dk' ± 9(x')9(l — x') f Amim2 a . 2 



x(l-x)x'(l -x') V 37r2 Q 2 



-X-W coni (Q 2 )jtfj(x',k' ± ) , (1) 

where M is the mass of the bound-state and ip is the projection of the light-front wave- 
function in the quark- ant iquark Fock-state. The confining interaction is included in the 
model by W con f(Q 2 ). The momentum transfer Q is the space-part of the four momentum 
transfer and the strength of the Coulomb-like potential is a. The singular interaction is 
active only in the pseudo-scalar meson channel with A as the bare coupling constant. 

For convenience the mass operator equation is transformed to the instant form represen- 
tation || , which in operatorial form is written as [f|] 

(M 2 + V + V 5 + V coni ) \<p) = M 2 \ip) , (2) 

where the free mass operator, M (= E\ + E 2 ), is the sum of the energies of quark 1 and 2 



(Ei = ymf + k 2 , i=l, 2 and k = \k\), V is the Coulomb-like potential, V s is the short-range 
singular interaction and Konf gives the quark confinement. 

We simplify Eq. by omitting the Coulomb term to the form |3| 

{Ml + g5(r)) cp(r) = M 2 tp(r) , (3) 
where the bare strength of the Dirac-delta interaction is g, and the mass squared operator 

is! 

ML = (C(k)k 2 + mf) + 2m s v(r) , (4) 
in units of H — c — 1, m s — mi + m 2 . The dimensionless factor of k is 

C (*) = 2 + f±^ + f±^. (5) 

In the following, we approximate C(k) as m s /m r 0, with m r = m^/^i + ^2)- 
The harmonic oscillator potential is introduced as a confinement 

v(r) = -c + ^c 2 r 2 , (6) 

where Co and c 2 are two universal parameters valid for all of the mesons. The eigenvalue 
Eq. (H) is given now by 

2m s (~^V 2 + l -c 2 r 2 + l -m s - c ) * n (f) = M 2 ^ n (f), (7) 

with ty n (r) the eigenstate of the harmonic oscillator potential and the corresponding eigen- 
value 

M 2 = 2m s (J^ n + \J c 2l m r + ~ c o) = nw + m 2 s + 3m s ^c 2 /m r - 2m s c , (8) 



where n (0, 1, 2, • ■ ■) is the radial quantum number and w = Am s Jc2/m r . Note that here 



n begins from for convenience, while in discussions in Section III, it begins from 1 to keep 
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accordance with literature. In the present model, Eq. @ gives the vector meson spectrum 
since Dirac-delta interaction acts only on pseudo-scalar mesons. The mass squared of the 
ground state (n = 0) of the 3 Si meson can be written as 



M: 



ml + 3m s ^C2/m r — 2m s CQ. 



(9) 



With the subtraction point, /i, taken to be the mass of the pseudo-scalar meson ground 
state, the reduced T matrix was derived in Ref. H as 



^ 1 (M 2 ) = (2 7 r) 3 El^(0)p 



1 



v /i 2 -Ml M 2 — Ml 
The value of the S-wave eigenfunction at the origin is given by [|1C 

"2 2 " n (2n + l)\ 
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where a 1 is the oscillator length. The final form the reduced T matrix is 

1 1 



nw — M? c 



M 2 — nw — M 2 S ( 



:i2i 



The zeros of Eq. (|T2|) gives the eigenvalues of the the mass squared operator of Eq. (|3|) . In 
Ref. ||, /i was changed continuously to study the splitting between the 1 Sq and 3 Si spectrum 
of the 7r-p mesons. In the following section, we are going to discuss this splitting in more 
detail and apply the model to mesons with heavy quarks as well. 



III. RESULTS AND DISCUSSIONS 
A. Parameters and nomenclature 

The masses of p(770), p(1450), J/i/j(1S) and i/j(2S) || and Eqs. (|, |9|) are used to fix c , 
C2 and up, down and charm quark masses with the assumption of m u = m^. Then masses of 
strange and bottom quarks are determined by the masses of K* and B* || and Eq. In 
order to predict t-quark meson spectrum, we adopt an estimate from Godfrey and Isgur , 
m t = 35 GeV. The parameters used in the present model are listed in Table |. 

The physical nomenclature of mesons are given in Table H to facilitate the following 
discussion. Since we assume that up and down quarks are the same, and since mixing 
between different flavors can not be dealt with within this simple model, we investigate only 
1 = 1 states among the diagonal meson sectors containing u, d and s quark. m u = 
also means that in our model the spectrum of vr° (p°, K°, K°\ D°, D *, B°, B *, T°, T *) 
is same as that of vr ± (p ± , i\T ± , i\T ± *, D ± , D±*, B ± , 5 ± *, T ± and T ± *). In the following 
discussions, we will omit the charge signs of mesons for simplicity. 
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B. S- wave meson spectra 



The splitting of the light mesons, ir-p and K-K* spectra, due to the Dirac-delta inter- 
action acting in the pseudo-scalar 1 5'o channel, are studied in the previous work ||. There 
the empirical slope (w) from Ref. || or Ref. || are used directly. In the present paper, ir-p 
and K-K* are reinvestigated in the same framework as other S'-state mesons. The results 
are presented in Fig. [I] and Fig. |2| respectively. Similar agreement with the data as that in 
Ref. [|| is found in Fig. [I], simply because the present model gives w = 1.55 GeV 2 which 
is very close to w = 1.39 GeV 2 used in Ref. [§]. For K*, the present model gives a large 
value of w, 1.92 GeV 2 , (due to the fact that strange quark mass is larger than up-down 
quark mass, w for K* must be larger than that for p from the present model), compared 
to that extracted from the data, 1.19 GeV 2 . However, as pointed out in Refs. there 
are ambiguities about the quantum number assignment of exited states of K and K*. If 
we follow the identification of the quark model § for #(1460) (2%), #(1830) (3%) and 
#*(1410) (2 3 «Si), the spectra of K-K* from present model are in good agreement with the 
data. 

In Fig. [3] the results are shown for the n c -if> mass splitting as a function of the ground-state 
pseudo-scalar mass p, which interpolates from the n c to the J / if) meson spectrum. Compared 
to the light ir-p and K-K* mesons, the splitting is smaller even for the ground state (~ 100 
MeV) and becomes weaker in the excited states, although the model attributes consistently 
smaller masses for the 1 S states compared to the respective 3 Si ones. An excited state of 
T) c is observed without definite spin and parity assignment |J. From our model, this state 
might be n c (2S), which is consistent with the assignment of quark model ||. There are 
many exited states for if>, such as ?/>(3770), ^(4040), ^(4160), and '0(4415), all of them are 
assigned to J n = 1". Considering that ^(3770) is ^(l 3 £>i) ^(4040) or ^(4160) seems 
to be ip(3 3 Si) from the present model. 

From Eqs. (|| |9[) and quark masses listed in Table |, one obtains the spectrum of D* . 
The splitting of D-D* and spectrum of D meson can be calculated from Eq. ([H]) with the 
renormalized strength of 5 potential fixed by the mass of corresponding pseudo-scalar ground 
state D and is shown in Fig. One finds good agreement for D*(1S) with the data. The 
predicted mass of D(2S) is about 10% larger than the unconfirmed data. 

Similarly, the present model gives D*(1S) with good agreement with the data as is seen 
from Fig. |^ in which the spectra for D s -D* s are presented. 

As mentioned before, the bottom quark mass is fixed by mass of B*(1S) and Eq. @. 
There are no more data for 6-quark mesons. The predicted spectra for B and B* are 
presented in Fig. || 

The spectrum of T can be calculated from our model with parameters listed in Table | 
and is shown in Fig. [7]. We should note that the agreement with the data is good, from 
T(IS') to T(4S"), considering that no parameter is adjusted specially for T, the heaviest 
meson states observed up to now. The unconfirmed bottomonium state 7^(9300) is used to 
predict the spectrum of rj b and also presented in the same figure. 

No confirmed data are available for the triplet states of B s and B c . With the bottom 
quark mass nib = 5068 MeV, the present model predicts the triplet ground state masses for 
B s and B c , being 5342.1 MeV and 6.346 GeV, respectively. Although the predicted mass for 
B S (1S) is comparable to the unconfirmed data 5416.6 MeV, both predictions are respectively 
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smaller than their corresponding singlet ground state masses, 5369.6 MeV and 6.4 GeV. 
Therefore we can not calculate within our model the singlet spectra B s and B c . However, 
for completeness and in order to give a reference to the reader, we list the calculated values 



for triplets of the bottom strange and bottom charged mesons in Table |TXT| and Table [TV 
respectively. 

The vector top meson spectra from the present model with the top quark mass m t = 35 
GeV M and m t = 175 GeV M are given in Table |V[ The top quark mass m t = 35 GeV 
is used in order to compare the present spectra with predictions of Ref. j7|. Qualitative 
agreement between the two models is found. 

C. Further discussions 

From the above results, it is clear that our model, with the mass squared operator 
consisting of a harmonic oscillator confining potential and a Dirac-delta interaction, could 
be used to describe universally and satisfactorily both singlet and triplet states of S'-wave 
mesons as well as radial excitations. However, lattice QCD calculations predict that the 
quark- ant iquark potential increases linearly with the distance between quark and antiquark 
r when r is large. One may ask, does this contradict our model? The answer is no. In 
the following, we justify roughly that the two interactions are consistent with each other for 
large r disregarding the case for small r which is not important in the present work. 

In the front form of QCD, the mass squared M 2 = m 2 + 2m s E = m 2 s + 2m s (TFF + Vff) 
(m s is the sum of the quark and antiquark mass, T is the kinetic energy and Vff is the 
interaction) ||, while in the instant form of QCD, M = m s + E = m s + (Tip + Vw). The 
relation between the interactions Vff and Vip is found for large r (where T < 7 < V 2 ) 
as Vff ~ V^p, from which the harmonic like potential for the mass squared operator in the 
front form can be derived from the linear confinement potential in the instant form of QCD. 

Very recently the radial excitations of light mesons were studied in detail in the framework 



of the QCD string approach [12]. There the spin-averaged meson masses were calculated 
with a modified confining potential and the calculated slopes of the radial Regge trajectories 
are in agreement with Ref. ||. In Ref. ||12|| , the linear relation between mass squared M 2 and 
the radial quantum number n comes mainly from properties of the approximated eigenvalues 
of the spinless Salpeter equation with the linear confining potential. Since in our light-cone 
QCD model, the harmonic oscillator potential is included in the mass squared operator, one 
arrives naturally at the same linear relation between M 2 and n for vector mesons. This 
gives an explicit explanation of the radial Regge trajectories found for light vector mesons 
in Ref. []5J . This relation is still valid even for heavy vector mesons as shown in the previous 
subsection. In addition, an extension of the present model with orbital excitation included 
could also be used to describe the orbital Regge trajectory. 

For pseudoscalar mesons, particularly for the light ones, the present model does not 
support the simple linear relation because the Dirac-delta interaction plays an important 
role now. However, this is not in contradiction with the data because light pseudoscalar 
mesons follow the radial Regge trajectories poorly. Let's take the pion as an example, the 
slope w = 1.67 GeV 2 is derived from the masses of it (M = 0.14 GeV) and its first radial 
excited state (M = 1.3 GeV), while a slope smaller by ~ 10%, w = 1.55 GeV 2 can be 
calculated from its first and second excited states (M = 1.8 GeV). 
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IV. SUMMARY 



We applied the renormalized light-cone QCD-inspired effective theory with confinement 
in the mass squared operator to study mesons with heavy quarks. For the confinement, 
the harmonic oscillator potential is used which allows an analytic solution of our model. 
The Coulomb-like potential is omitted in the present work while the Dirac-delta interaction 
is kept which acts on the singlet S'-wave states and plays important role in the splitting 
of the singlet-triplet spectra. The two parameters of the harmonic potential and quark 
masses are fixed by masses of p(770), p(1450), J/ip, if)(2S), K*(892) and B*. A T matrix 
renormalization method is used to renormalize the model, in which the pseudo-scalar ground 
state mass fixes the renormalized strength of the Dirac-delta interaction. 

The model is applied to study the 5-wave meson spectra from n-p to rjb-T and is also used 
to predict top quark meson spectra. The linear relationship between the mass squared of 
excited states with radial quantum number — the radial Regge trajectory [|T2J — is apparent 
from our model and is found to be qualitatively valid even for heavy mesons like T. 

The simple model presents satisfactory agreement with available data and/or with the 
meson mass spectra given by Godfrey and Isgur 0. Therefore, the recently proposed ex- 
tension of the light-cone QCD-inspired model which includes confinement while keeping 
simplicity and renormalizability, gives a reasonable picture of the spectrum of both light 
and heavy mesons. An extension of the present model with orbital excitations included 
could also be used to describe the orbital Regge trajectories. 
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TABLES 



TABLE I. Parameters used in the present paper, cq and C2 of harmonic oscillator potential 
and masses of up, down and charm quarks are fixed from masses of p(770), p(1450), J/tp(lS) and 
tp(2S), with the assumption of m u = m^. Strange and bottom quark masses are determined by 
masses of K* and B*. The top quark masses, rat = 35 GeV, same as that in Godfrey et al. 0, 
and nit = 175 GeV from Q, are used to predict spectra for i-quark mesons. The data for meson 
masses are taken from Hagiwara et al. ||. 

Parameters c [MeV] c 2 [GeV 3 ] m u = m d [MeV] m s [MeV] m c [MeV] m b [MeV] m t [GeV] 
Values 807 7.13xl0~ 2 265 478 1749 5068 35/175 



TABLE II. The physical nomenclature of the mesons as a reminder. Pseudo-scalar mesons are 
given on the left, vector mesons on the right of each sector. Diagonal sectors are marked with 
frames for guiding eye. 



d 



'•'/•'/' p" 



7T 



7i ||p' 

K°\\K*° 
D + \\D* 

I?\\B* 

rp-\- I I rp* 



k°\\k*° 

K + \\K*+ 



K-\\K*- 
D°\\D*° 
B~\\B*- 

^0 1 1^*0 



D-\\D*- 
D°\\D*° 
D7\\D*~ 



Df\\D*+ 

R°llR*° 
B s\\ B s 

t+\\t: + 



B°\\B*° 
B + \\B* + 
B° s \\Bf 
B+\\B*+ 



Bc\\B* c - 



m. 



T 



± c I \ ± c 



*0 



1 b W 1 b 



7^0 

c 



rp* — 

rp* — 
S 

c 



TABLE III. The present spectra of vector bottom strange mesons B*(nS) compared with 
available data from Hagiwara et al. || and predictions of Godfrey et al. |?J . Note that the data for 
B*(1S) is not confirmed J7j. Masses are in MeV. 



B* Ours Data §] Godfrey et al. § 



l 3 Si 5342.2 5416.6 5450 

2 3 Si 6123.8 6010 

3 3 Si 6816.4 

4 3 Si 7444.7 

5 3 Si 8024.2 

6 3 Si 8564.6 
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TABLE IV. The present spectra of vector bottom charmed mesons B*{nS) compared with 
predictions of Godfrey et al. |7j (no data available). Masses are in MeV. 



Bl Ours Godfrey et al. § 



2 3 S 1 6830.4 6890 

7282.9 



PSi 6345.8 6340 

4 3 Si 7709.0 

5 3 Si 8112.6 

6 3 Si 8497.1 



TABLE V. The present spectra of vector t-quark mesons compared with predictions (in paren- 
thesis) of Godfrey et al. |?J. Results with m t = 35 GeV [|7| and m t = 175 GeV || are given in the 
first and the third row of each sector. The masses are in GeV. The spectra of T*° and T*^ are the 
same since m u = ntd is assumed. 



Mesons 


M(l 3 5i) 


M(2 3 Si) 


M(3 3 5i) 


M(4 3 Si) 


M(5 3 5i) 


M(6 3 Si) 


rp* 


35.24 


36.27 


37.27 


38.24 


39.19 


40.12 




(35.39) 
175.24 


176.27 


177.30 


178.33 


179.34 


180.36 


1 S 


35.25 


36.03 


36.79 


37.53 


38.26 


38.97 




(35.46) 
175.25 


176.02 


176.79 


177.56 


178.32 


179.08 


1 C 


36.25 


36.67 


37.08 


37.49 


37.89 


38.29 




(36.30) 
176.25 


176.65 


177.06 


177.46 


177.87 


178.27 


n 


39.45 


39.71 


39.96 


40.21 


40.47 


40.72 




(39.31) 
179.44 


179.68 


179.92 


180.16 


180.40 


180.64 


9 


69.29 


69.42 


69.54 


69.67 


69.80 


69.93 




(68.70) 
349.24 


(69.39) 
349.29 


(69.71) 
349.35 


(69.93) 
349.41 


(70.10) 
349.46 


(70.40) 
349.52 
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FIG. 1. Mass of the excited qq states (M*) as a function of the mass (//) of the pseudoscalar 
meson ground state for 1 = 1. The data is taken from Hagiwara et al. Q and is shown by solid 
circles on the left and right of the figure for 1 5'o and 3 S\ mesons, respectively. p(1700) (labelled 
with empty diamond) might be D-wave dominant PJTl||. Calculated 1 So and 3 Si spectra are given 
in parentheses within the dashed lines. 
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FIG. 2. Mass of the excited qq states (M*) as a function of the mass (p) of the pseudoscalar 
meson ground state for / = 1/2 of the strange mesons. The data is taken from Hagiwara et al. || 
and is shown by solid circles on the left and right of the figure for 1 So and 3 Si mesons, respectively. 
Errors are not available for iT(i460) and i^(1830). There is ambiguity about the radial quantum 
number of if*(1410) and iT*(1680) || (labelled with empty diamond). -fT(3100) is not confirmed 
and represented by empty square. Calculated 1 5o and 3 Si spectra are given in parentheses within 
the dashed lines. 
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FIG. 3. Mass of the excited qq states (M*) as a function of the mass (//) of the pseudoscalar 
meson ground state for I = of the charmed mesons. The data is taken from Hagiwara et al. || 
and is shown by solid circles on the left and right of the figure for 1 So and 3 Si mesons, respectively. 
No radial quantum numbers assigned to -0(3770), ^(4040), ^(4160) and ^(4415) (labelled with 
empty diamond). r/ c (3594) is not confirmed and represented by empty square. Calculated 1 So and 
Si spectra are given in parentheses within the dashed lines. 
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FIG. 4. Mass of the excited qq states (M*) as a function of the mass (//) of the pseudoscalar 
meson ground state for / = 1/2 of the charmed mesons. The data is taken from Hagiwara et al. || 
and is shown by solid circles on the left and right of the figure for 1 Sq and 3 Si mesons, respectively. 
.D*(2637) is not confirmed and represented by empty square. Calculated 1 S'o and 3 Si spectra are 
given in parentheses within the dashed lines. 
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FIG. 5. Mass of the excited qq states (M*) as a function of the mass (p) of the pseudoscalar 
meson ground state for / = of the charmed strange mesons. The data is taken from Hagiwara 
et al. H and is shown by solid circles on the left and right of the figure for 1 So and 3 Si mesons, 
respectively. Calculated 1 S'o and 3 Si spectra are given in parentheses within the dashed lines. 
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FIG. 6. Mass of the excited qq states (M*) as a function of the mass (p) of the pseudoscalar 
meson ground state for / = 1/2 of the bottom mesons. The data is taken from Hagiwara et al. || 
and is shown by solid circles on the left and right of the figure for 1 Sq an d 3 S± mesons, respectively. 
Calculated 1 S'o and 3 S± spectra are given in parentheses within the dashed lines. 
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FIG. 7. Mass of the excited qq states (M*) as a function of the mass (p) of the pseudoscalar 
meson ground state for I = of the bottom mesons. The data is taken from Hagiwara et al. || 
and is shown by solid circles on the left and right of the figure for 1 Sq and 3 Si mesons, respectively. 
No radial quantum numbers assigned to Y( 10865), and T(11020) (labelled with empty diamond). 
7/5(9300) is not confirmed and represented by empty square. Calculated x 5o and 3 Si spectra are 
given in parentheses within the dashed lines. 
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